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Dynamic analysis of aerospace systems is required to in- 
sure that the structures will maintain their integrity 
and provide predictable and acceptable mechanical per- 
formance throughout the mission profile of specified 
acceleration environments. This analysis includes the 
structural response to shock and vibration and evaluates 
the maximum deflections and material stresses and the po- 
tential for the occurrence of elastic instability, fa- 
tigue and fracture. The required computations are often 
performed by means of finite element analysis (FEA) com- 
puter programs in which the structure is simulated by a 
finite element model which may contain thousands of ele- 
ments. The formulation of a finite element model can be 
time consuming, and substantial additional modeling ef- 
fort may be necessary if the structure requires signifi- 
cant changes after initial analysis. This paper pre- 
sents rapid methods for obtaining rough estimates of the 
structural response to shock and vibration, for the pur- 
pose of providing guidance during the intial mechanical 
design configuration stage. 


INTRODUCTION 

Structures are often made up of simple components such as beams, rings, 
arcs, plates and shells. The natural frequencies of such a structure cannot 
usually be found from the frequencies of these components. However, the stiff- 
ness, damping and mass of these components, the stiffness and damping of the 
connections between components, and the type of attachment of the structures to 
mounting surfaces will determine the natural frequencies of the structure. 
Estimates of natural frequencies can only be made for simple structures without 
developing an FEA (finite element analysis) model and utilizing an FEA computer 
program. But even rough estimates of natural frequency can provide a rela- 
tively rapid means of comparing maximum acceleration, stress and fatigue in dif- 
ferent design approaches, and identifying potential problem areas in a struc- 
ture. This type of information can help to avoid excessive modification of the 
FEA model when a comprehensive computer analysis is done. 

The methods presented here for evaluating natural frequencies of simple 
structures will typically have a frequency error on the low side, which will 
result in a conservative (larger than actual) estimate of stress. 
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Composite Beams 

Beams may be made of two or more layers of different materials adhered to 
one another, with each layer running the length of the beam. In this section, 
the layers are assumed to have constant, rectangular cross sections. The 
layers may be oriented so that the direction of vibration is parallel to the 
layer interfaces or normal to the layer interfaces. 


Vibration Parallel to Layer Interfaces 


Figure 1 shows the case where the direction of vibration is parallel to 
the layer interfaces. The term EI/L is referred to as the stiffness of a beam. 
Defining El as the stiffness factor of a beam, the stiffness factor of the com- 
posite beam of Figure 1 is 


n n 

El - Z Ejli = (1/12) £ Eibihi 3 lbf.in 2 
i=l i=l 


( 1 ) 


where E^ 
bi 
hi 
li 


modulus of elasticity of layer i, lbf/in 2 , 
width of layer i, inch, 
height of layer i, inch, 

area moment of inertia of layer i about neutral (Z) axis, in • 
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Fig. 1 Layered Beam, Vibration Parallel To Layer Interfaces 
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The weight per unit length of the composite beam is 


" ( 2 ) 

W = 2 p i b-j h-j Ibf/in 
i=l 

where P, = weight density of layer i, lbf/in^ 

Equations (1) and (2) may be used with slender beam frequency formulas when the 
composite beam is uniformly loaded and the neutral axis through the beam cross 
section, parallel to the Z axis, remains undeflected (no bending along the Z 
axis). 

Vibration Normal to Layer Interface 

Figure 2 shows the case where the vibration is normal to the layer inter- 
faces. The stiffness factor of the composite beam of Figure 2 is 


El = Z b i h i E i [(Y-Y i )2 + ( l/12)h-j2] lbf.in2 
i=l 

i-1 

where Y-j = 2 hj + ( 1/2 ) h-j inch, where h 0 = 0, 
j=0 

__ n / n • 

and Y= 2 b-j h-j E-j Y n - / 2 b-jh-j E-j inch 
i=l / i=l 

See Equation 1 for definitions. 



Fig. 2 Layered Beam, Vibration Normal to Layer Interfaces 
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The weight per unit length of the composite beam, W, is given by Equation 
Equations (2) and (3) may be used with slender beam frequency formulas to 
obtain approximate values of natural frequency. Accuracy is improved when the 
layer widths, b^ , approach equality with one another. 


Stepped Beams 


Stepped beams have two or more different cross sections along their span, 
resulting in two or more different moments of inertia. Figure 3 shows two 
examples of stepped cantilever beams. In Figure 3a the beam has two different 
cross sections and the average moment of inertia for the beam is 


I A = L 3 l 1 I 2 /[3(a 2 b+b 2 a+b 3 /3)I 1 + a 3 I 2 ] in 4 ( 4 ) 

In Figure 3b the beam has three different cross sections and the average moment 
of inertia for the beam is 


I A = L 3 I 1 I 2 l 3 ^ <^3[(a+b) 2 c + (a+b)c 2 + c 3 /3]I 1 I 2 
+ 3(a 2 b+ab 2 +b 3 /3)I^l3 + a 3 I 2 l 3 ^ in 4 


(5) 
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Fig. 3 Stepped Cantilever Beams With a) Two Different Cross Sections, and 
b) Three Different Cross Sections 
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In general, for a beam with n different cross sections, an approximate value 
for the average moment of inertia is 


n 

1^. = (1/L) 2 x^I^ in^ 
i = l 

where x A = spanwise length of cross section i, inch. 

Ii = moment of inertia of cross section i, in^, 
L = 2 x£ = full span of beam, inch, 
i 


( 6 ) 


Equations (4), (5) and (6) may be used in standard frequency formulas to obtain 
values of natural flexural frequencies of stepped beams. Equations (4) and (5) 
are only for cantilever beams and should provide accurate results for slender 
beams. Equation (6) may be used for any end support conditions, and will usu- 
ally yield a natural frequency roughly 5% to 10% lower than the correct value. 


Slender Right Angles and U Bends 


Figure 4 shows a right angle and a U bend with intermediate supports. The 
ends, E, may have any combination of pinned, P, or clamped, C, boundary con- 
ditions. The intermediate supports, S, prevent transverse motion (perpendic- 
u ^ ar the beam axis) at the support but allow the beam to move parallel to 

its own axis and to rotate about any axis. The fundamental natural frequency 
for vibration in the plane of the figures (in-plane vibration) is [2] 

F = (A/2TTR2) • (EI y g/W) ^ Hz (7) 

where A = dimensionless frequency parameter in Table 1, 

R = radius of curvature shown in Figure 4, inches, 

E = modulus of elasticity of beam material, lbf/in 2 , 

Iy = area moment of inertia about axis perpendicular to the plane of 
the figures, in^, 

g = gravitational acceleration at surface of earth = 386 in/sec 2 , 

W = weight per unit length of beam, lbf/in. 


The fundamental natural frequency for vibration perpendicular to the plane of 
the figures (out-of-plane vibration) is 

F = (X/27TR 2 ) • (Glpg/W)^ Hz (8) 


where R, E, g and W are as defined for Equation (7), 


A - dimensionless frequency parameter in Table 2, 

G = shear modulus = E/2(l+v) lbf/in^, 

V = Poisson’s ratio, dimensionless, 

Ip = polar area moment of inertia = I x +I y , in 4 , 

I x = area moment of inertia about axis in the plane of the figure 
and perpendicular to the local beam axis, in 4 , 

Iy is as defined for Equation 7. 
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Fig. 4 Slender a) Right Angle Bend, and b) U Bend 
Table 1 In-Plane Vibration of Right Angles and U Bends 



NOTES: 1. L = length of legs, inch 

2. R = radius of curvature, inch 

3. P = pinned end condition 

4. C = clamped end condition 

5. Vibration in the plane of Figures 4a and 4b 


I 

i 
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Table 2 Out-of-Plane Vibration of Right Angles and U Bends 



NOTES: 1. See notes 1-4 of Table 1 

2. Vibration perpendicular to the plane of Figures 4a and 4b 

3. x values for v = 0.3 


Equations (7) and (8) do not take into account shear deformation, cross 
sectional distortion due to torsion, or coupling of rotation and displacement. 
The rotary inertia of the beam twisting about its own axis is not included in 
Equation (7) but is included in Equation (8). However, the values of A. given 
in Table 2 are only valid for circular beams or tubes with a value of V =0.3. 

Simple Frames 

The simple frames shown in Figure 5 are also called portal frames in struc- 
tural applications or bents in electronic applications. The following formulas 
provide approximate values for the fundamental natural frequencies in the 
specified vibration modes. 



DIRECTION 

OF 



Fig. 5 A Simple Frame in a) In-Plane Vertical Vibration, 
b) In-Plane Lateral Vibration, and c) Out-of-Plane Transverse Vibration 
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5a, 


For in-plane vertical vibration with legs hinged at the supports, Figure 


F = (1/2TT) ^A8EI 1 g/WL 3 [l-2.25/(2K+3)]J' 1 Hz 

where E = modulus of elasticity of frame material, lbf/in 2 , 

1^ = area moment of inertia of top of frame about neutral axis, in , 
g = gravitational acceleration at surface of earth = 386 in/sec , 

W = total weight of frame, lbf , 

L = length of top (span of frame), inch, 

K = hIj/Ll 2 , dimensionless, 
h = height of frame (length of leg), inch, 

12 = area moment of inertia of frame legs about neutral axis, in . 


(9) 


5a, 


For in-plane vertical vibration with legs fixed at the supports. Figure 


F = (1/2tt) |48EI 1 g/WL 3 [l-3/(2K+4)]J ** Hz 
For in-plane lateral vibration with legs 
F = (1/20 j^24EI 2 g/Wh 3 [l+3/(6K+l)]j^ Hz 


( 10 ) 


For in-plane lateral vibration with legs fixed at the supports, Figure 5b, 

(ID 


For out-of-plane transverse vibration with legs fixed at the supports, 
Figure 5c , 


F = (g^/2u) < (W/2) [L 3 /24EI 1 + h 3 /3EI 2 
- L 4 GC 2 /32EI 1 (2hEl 1 +LGC 2 )] > Hz 

where C 2 = torsional constant, inO 


( 12 ) 


The approximate fundamental natural frequency for a rigid body of mass M Q 
supported by n slender, uniform legs of length L, all in the same plane, clamp- 
ed at their feet and at the rigid body, as shown in Figure 6, for vibration in 
the plane of the legs, is given by^l 

F = (1/270 [(12^E i I i )/L 3 (M o +0.372M i )] * Hz (13) 


where M 0 = rigid body mass, lbs. mass, 

= mass of leg i, lbs. mass, 

Ei = modulus of elasticity of leg i, lbf/in 2 , 

=> area moment of inertia of leg i about its neutral axis, in , 

E = sum over all legs, i = 1, 2, 3, ..., n, 
n = number of legs 2. 

Housings 

Housings may be analyzed to estimate their fundamental natural frequen- 
cies. These frequencies may include flexural vibration along one or more axes 
of the structure, torsional vibration, and coupled modes of vibration. The 
frequencies will depend on the geometry and material properties of the 


210 



DIRECTION OF 
. VIBRTATION . 


M o y 

M 1 

m 2 

m 3 

• • 

El 

e 2 

E 3 


h 

>2 

1 3 

• • 

■ — ... • • 



7777 - 7777 - 7777 


Fig. 6 Rigid Body on Slender Legs 

structure f the attachment efficiency factor between parts of the structure 

and the connection of the structure to mounting surfaces (the boundary condi- 
tions). 

Flexure 


Figure 7 shows a housing composed of several structural elements. It is 
mounted by means of brackets attached at each end of the longer dimension, near 
the bottom (Figure 7b). It will be analyzed as a simply supported beam which 
can vibrate in flexure in the X direction and in the Y direction, and which can 
vibrate in torsion about the Z axis. Evaluation of the flexural frequencies 
requires an estimate of the stiffness factors of the structure, E I and El. 
These are determined as follows. x x Y Y* 


6 _ 

E X I X = ^ riiEiUitX-Xi) 2 + I x>i ] 
i=l 


lbf 


(14) 


where 


n i = attachment efficiency factor for element i dimensionless, 

E i = modulus of elasticity of element i, lbf/in*, 

A i = cross sectional area of element i in the x-y plane, in 2 , 

x i ~ distance from left edge of structure to neutral axis (or 

mid-point) of element i, inch, 

Ix,i “ area moment of inertia of element i about the neutral 
axis parallel to the Y direction at X ± , in 4 , 
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B) FRONT VIEW 


Fig. 7 


2 



and 


__ 6 6 

X = £ E ri^A^E^ inch, 

i=l i=l 


(15) 


E y x y 


6 

Z ^EitAiCY-Yi) 2 + I t ] lbf.in 2 
i = l 7 * 


(16) 


where nj.Ej and are defined following Equation (14), 


Y i - distance from bottom of structure to neutral axis (or 
mid-point) of element i, inch, 

x y,i = ar ea moment of inertia of element i about the neutral axis 
parallel to the X axis at y it in 4 , 


and 

6 6 

Y * ^ ^i^i^i^i/ Z ^i^i^i inch 

i“l i=l 


(17) 


Note that structural elements 7, 8 and 9 are not included in Equations (14) 
through (17) because they are not subjected to bending but are either fixed to 
the mounting brackets or displaced parallel to their own plane. The funda- 
mental natural frequencies for flexural vibration in the X and Y directions are 

F x = (TT/2L 2 )(E x I x g/W) 15 Hz (18) 

and 

Fy = (Tr/2L 2 )(EyIyg/W) Ji HZ 

where L = length of housing, Figure 7b, inches, 

g * acceleration of gravity at surface of earth = 386 in/sec 2 , 

W = weight of housing per unit length, lbf/inch (Total weight = WL) 

E X I X and E y I y are found from Equations (14) and (16). 

If the ends of the housing, structural elements 7 and 8, were mated with and 
fixed to mounting surfaces, then the constant in Equation (18) would be 
(22.373/tt) instead of it. 

Torsion 

Acceleration in the X direction will produce torsion as well as bending 
since the housing is supported near its bottom and the center of gravity (CG) 
is located above the support (Figure 7b). The axis of rotation for torsion 
will be at the bottom of the housing parallel to the Z axis. See Figure 8. 
Since the CG is not on the axis of rotation, the torsional natural frequency of 
the housing will be coupled to the flexural natural frequency, F x . 
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Fig. 8 Torsional Deflection, With One-Half of Housing Shown 


The frequency formula for torsional vibration is 

F 0 = (A/2TTL) • (CGg/yl p ) ^Hz (19) 


where 


C = torsional constant of beam cross section, in^ 

G = shear modulus of beam material, lbf/in^ 

U = weight density of beam material, lbf/in 3 

Ip = polar area moment of inertia of beam cross section about the beam 
axis of torsion, in^. 


The polar mass moment of inertia, J, about the axis of rotation is an important 
factor in the present case. The half-housing shown in Figure 8 may be analyzed 
with the frequency constant X found by solving the transcendental equation 
which includes the parameter J: 

cotX = (Jg/yLC)X (20) 

For the half-housing, L must be replaced by L/2 in Equations (19) and (20), and 
J must be the polar mass moment of inertia about the axis of rotation for the 
half-housing. The same results may be obtained by using the formula 

F q = (1/2IT)* (4CG/LJ )** Hz (21) 

U 

instead of Equation (19) when (ULIp/gJ) 2 << 1. In Equation (21), L is the 
full length of the housing and J is the polar mass moment of inertia about the 
axis of rotation for the full housing. C and G have the same values as in Equa- 
tions (19) and (20). Using the relation U = Mg/LA lbf/in 3 , where M is the 
total mass of the housing and A is the housing cross sectional area in the X-Y 
plane, and equating Equations (19) and (21) yields: 
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Ip/A = X 2j/ M in 2 (22) 

where the radius of gyration is (Ip/A)^. In arriving at Equation (22), L was 
replaced by L/2 in Equation (19). 

The value of the polar area moment of inertia about the axis of rotation. 
Ip, may be estimated by 

*p = (E X I X + Eyly)/E a in* (23) 

where E a is an appropriate average value for the structure. In the case where 
the structure is composed primarily of a single structural material, E a is the 
modulus of elasticity for that material. 

The value of J is given by 


J - £ [Jz,i + m i( x i^ + yi^)] lbmass.in^ 

i=l 


where J 2>i = polar mass moment of inertia about axis through the neutral 
axis of element i and parallel to the axis of rotation, 
lbmass.in^, 

m^ = mass of element i, lbmass, 
x i»yi previously defined. 

A more rapid but less accurate approximation for J is 

J = (M/12)(4h2+d2) lbmass. in^ (25) 

where M is the total mass of the housing. (See Figure 7.) 

If J is calculated for the entire housing, then only half of its value must be 
used in Equation (20). 


The torsional constant C is much more difficult to accurately estimate, 
even for a simple structure. The following rough approximation may be used:!^] 

C s [l/2]I p in* (26) 

Coupled Modes 

As pointed out previously, acceleration in the X direction produces both 
bending at a frequency F x and torsion at a frequency F q. These vibration modes 
will be coupled to produce a fundamental natural mode of the structure which 
can be approximated by Dunkerley's method 17] 

F c = (F x -2 + F 0 -2)- 1 ' 2 Hz (27) 

where F x and Fe are found from Equations (18) and (19). The coupling of modes 
of vibration always results in a natural frequency lower than the coupled 
frequencies, and a consequent increase in deflection and stress. In the 
example given in the preceeding two sections, this coupling of modes may be 
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avoided by mounting the housing so that the center of gravity (CG) lies on the 
mounting plane. Torsional modes may still occur but they will not be coup- 
led with the bending mode. Figure 9 is an example of a CG mount. 



Fig. 9 A CG Mount 


Other Housing Configurations 

There are other housing geometries and mounting configurations where the 
housing may be modeled as a beam or a plate with boundary conditions which ap- 
proximate the mounting attachments. In these cases the standard frequency form- 
ulas may be used to estimate the fundamental natural frequency of the housing. 
The flexural stiffness factors, El, and the torsional frequency parameters C, 
Ip and J must be estimated as in Section 5.3. 

The examples shown in Figure 10 could be analyzed as follows. 

Figure 10a. 

1) For L > h,t : A beam with simply supported ends. 

Flexural vibration in X and Z directions. 

Torsional vibration about the axis through centroid, parallel to Y 
axis . 

Flexural and torsional modes not coupled. 

2) For L s h > t: A plate with two opposite sides simply supported, other 
two sides free. 

Flexural vibration in Z direction. 
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COORDINATES: 


COORDINATES: 













d) 


h“ L — H 





H = HOUSING 
P = MOUNTING PLATE 
B = BRACKET 




Fill lit] 


3) For L = t > h: A plate with four corner supports. 

Flexural vibration in X direction. 

Figure 10b. 

For L s h > t: A plate with three simply supported sides. 

Flexural vibration in Z direction. 

Figure 10c* Same as 10b. 

Figure lOd. 

1) For h > L,t: A cantilever beam. 

Flexural vibration in Y and Z directions. 

Torsional vibration about the axis through centroid, parallel to X 
axis . 

Flexural and torsional modes coupled if CG is not at the mid - point of 
the L and the t dimensions. 

2) For L = h > t: A plate simply supported on one side with other three 
sides free. 

Flexural vibration in Z direction. 

Torsional vibration about the axis through centroid parallel to X axis. 
Flexural and torsional modes coupled if CG is not at the mid-point of 
the L dimension. 

Lumped Elements 

One type of model that may be used to represent structures is the lumped 
element model. In this approach, parts of the structures are treated as masses 
and other parts as springs. Rigid, heavy components may be treated as masses, 
while flexible, light weight components may be treated as springs. The spring 
elements may be some of the same structural components that make up the masses, 
even though they are treated as massless in the analysis. The combined weight 
of the masses must add up to the total weight of the complete structure. 

Figure 11 shows a three degree-of-f reedom structure composed of a trans- 
former mounted on a bracket which is attached to a PWB mounted in a housing. 
The bottom of the housing is fixed to a mounting plate. When the acceleration 
is parallel to the mounting plate and normal to the PWB, the housing will vi- 
brate as a cantilever, the PWB will vibrate as a loaded plate, and the brack- 
et will vibrate in the direction shown. Figure 12 is a lumped element model of 
the structure shown in Figure 11. The values of the spring constants, K, may 
be determined by use of the following equation: 

K = (w/g) • (2TTF) 2 lbf /inch (28) 

where W = weight of element, lbf, 

g = acceleration of gravity at surface of earth = 386 in/sec^, 

F = fundamental natural frequency of the element, Hz. 

The frequency formulas for finding the values of F in Equation (28) are given 
by Equation (9) for the bracket and by standard formulas for the PWB and the 
cantilever housing. The static deflections, for a one-g acceleration, are 
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Fig. 11 Example of Three Degree-of-Freedom Structure 
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Fig. 12 Lumped Element Model of Three Degree-of-Freedom Structure 


219 






( 29 ) 


Yl - (W]^+W 2 + W 3 )/K^ inch 
y 2 = y^+(W 2 +W 3 )/K 2 inch 


Y 3 = y^+y 2 +W3/K3 inch 

The fundamental natural frequency of the structure shown in Figures 11 and 12 
is t ® ^ 

3 3 ^ (30) 

F = (l/2ir)[g Z ViYi/ £ W i y i 2 ]' 5 Hz 
i=l i=l 

For a structure with n degrees of freedom, which can be represented by n 
spring/mass elements in series, 

i-l n (3D 

yi = Z yj + Z Wj/Kj inch 
j=l j=i 


and 


F = ( 1 / 2tt ) [ g £ W iyi / £ W iyi 2 ] ^ Hz 
i=l i-l 

Accelerations and Stresses 

The natural frequency of a structure will determine what will be its accel- 
eration response in a shock or vibration environment. Two different structural 
designs with two different natural frequencies will generally experience differ- 
ent accelerations, resulting in different stresses when exposed to the same 
shock or random vibration spectrum. The estimates of natural frequency allow 
comparison of structural designs in terms of dynamic loads due to the specified 
acceleration environments. 

Figure 13 shows the shock acceleration response spectrum for an NSI ordi- 
nance, one-inch separation nut. The acceleration response is the three-sigma 
peak value (exceeded only 0.28% of the time) in units of the gravitational 
acceleration at the earth's surface, and it peaks at 15,000g and 10,000Hz. A 
structure with a natural frequency of 200Hz will experience a 215g shock, while 
one with a natural frequency of 250Hz will experience a 310g shock. The 
acceleration ratio is 310/215 = 1.44. 

Figure 14 shows a typical random vibration power spectral density (PSD) to 
which a space structure would be exposed during launch. The three-sigma peak 
acceleration response is approximated by 

G pk = 3[(ir/2) • PSD • Q 0 * F 0 ] ^ g (33) 

where Q 0 = transmissibility at the structure's natural frequency, F Q . Assuming 
the same value of Q 0 for the two structures with natural frequencies of 200Hz 
and 250Hz , the acceleration ratio will be 1.12. 
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FREQUENCY, Hz 


Fig. 14 A Typical Random Vibration Power Spectral Density (PSD) 



In the above example, it is seen that the accelerations, and consequently 
the dynamic loads, between different structural designs may be compared when 
estimates of natural frequency can be made. 
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